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§ '■ ABSTRACT 

(N : 

! We present spectra of the extended ^^CO and ^^CO J =1^0 emission along 

^ ! the major axes of 17 nearby galaxies. Spatial variations in the ratio of CO and 

•/^ ! ^^CO integrated intensities, TZ, are found in nearly every galaxy observed. There 

^ i is an overall variation in TZ of 20-40% from the inner 2 kpc to the disk. Roughly 

^ ! a third of the survey galaxies have such gradients in TZ detected above the 2a 

i confidence level. Though some galaxies show a lower central value of TZ, on 

^ i average TZ inside 2 kpc is 10-30% higher than TZ outside of 2 kpc. The average 

O ! CO/^^CO intensity ratio within the central 2 kpc of the survey sources is 11.6±0.4 

^ I (based on the noise) ±1.5 (based on systematic uncertainties estimated from 

"r^ I daily variations in CO and ^^CO intensities). The la dispersion in TZ between 

Ph| galactic nuclei of 4.2 is also quite large. The average value of TZ outside 2 kpc is 

O ' 9.8 ± 0.6 ± 1.2 with a standard deviation of 4.5. 

I An increase in the CO/^^CO intensity ratio from disk to nucleus may imply 

I that the conversion factor between CO intensity and H2 column density, X, is 

.5^ I lower in galactic nuclei. Also variations in physical conditions, most notably the 

^ ' gas kinetic temperature, affect both TZ and X. Abundance variations probably 

_c3 I do not cause the gradient in TZ, though we do not rule out a decrease in effective 
cloud column densities in galactic nuclei possibly caused by destructive starburst 
superwinds. A modest rise in temperature (less than a factor of 2 or 3) from 
outside a 2 kpc radius towards the nucleus can easily account for the observed 
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gradient. These results support previous work implying that X is lower in the 
center of the Milky Way and probably most galactic nuclei. Therefore calculating 
H2 masses using the standard Galactic X-factor, especially within the central few 
kpc of galaxies, overestimates the true mass by factors of a few. The standard 
X-factor still appears to be appropriate for galactic disks. 

Subject headings: galaxies: ISM, starburst — ISM: clouds, molecules 

1. Introduction 

Central to a complete picture of galaxy evolution is understanding star formation over 
large scales. Given that stars form in molecular gas, the evolution of a galaxy depends on its 
gas distribution. The main constituent of this gas is molecular hydrogen, which, due to its 
symmetry, has no dipole transitions between low-lying levels (i.e., no permitted rotational 
transitions) . Consequently, the rotational lines of carbon monoxide and its isotopic variants 
(e.g., ^^CO and C^^O) have been used to study molecular gas in galaxies. The miUimeter- 
wave lines of CO are the brightest among non-masing molecular transitions. In particular, the 
J=l— i^O line of ^^CO (hereafter written simply as CO) has been invaluable for determining 
the distribution of molecular gas column densities within the Milky Way and other galaxies 
(Dame, Hartmann & Thaddeus 2001; Young et al. 1995). 

These distributions are inferred assuming that the molecular gas column density, Nu^, 
is traced by the velocity-integrated radiation temperature (Rayleigh- Jeans brightness tem- 
perature), /(CO J = l— >0) = jTjidv. Though the CO emission is usually optically thick 
and consequently may not reliably trace column densities, the ratio Nn^/Ico, oi "X-factor," 
is found to be roughly constant on the size scales of many parsccs, with a value of X = 
1.6x10^° cm~^ (K kms^^)^^ in the Galactic disk (Dame, Hartmann & Thaddeus 2001; 
Hunter et al. 1997). The conventional explanation of how the optically thick CO line consis- 
tently measures column density requires that giant molecular clouds are in self-gravitational 
equilibrium and that the mean density and temperature do not significantly vary from cloud 
to cloud (Dickman, Snell & Schloerb 1986; Sakamoto 1996). Despite these results, there 
is evidence that X is 2-10 times lower than the standard value in the central few hundred 
parsecs of the Galaxy (Sodroski et al. 1995; Dahmen et al. 1998). It may also increase 
by 2-4 times the standard value in the outer Galaxy (Digel et al. 1990; Sodroski 1991; cf.. 
Carpenter et al. 1990). Nonetheless, a secondary measure of gas column density is required 
to gauge any variation of X with galactocentric radius. 

There are suggestions in the literature both for and against the X-factor varying with 
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position in external galaxies, often from observations of the rotational lines of ^^CO, a rarer 
isotopic variant of the more commonly observed CO. Since ^^CO is 30-70 times less abundant 
than CO (e.g., Langer & Penzias 1990), its rotational lines are usually optically thin on 
parsec scales. Therefore, assuming similar molecular gas physical conditions, observations of 
^^CO can give reliable molecular gas column densities, and any variation in the ratio of CO 
and ^^CO J=l— >0 integrated intensities, TZ, can then test whether X varies systematically 
within galaxies. Rickard & Blitz (1985) claim that this ratio is up to a factor of 5 higher in the 
central 1' than in the disks of six galaxies, implying that X is up to 5 times smaller in galaxy 
centers. In contrast, the CO/^^CO intensity ratios observed by Young & Sanders (1986) and 
Sage & Isbell (1991) show no significant differences between the centers of spiral galaxies 
and their disks. However, any study of molecular gas column densities that uses ^^CO lines 
must address the uncertainty in assumptions such as the ^^CO abundance, the gas kinetic 
temperature, and whether local thermodynamic equilibrium (LTE) applies. Also, 7^ is a 
sensitive function of optical depth. Hence, any observed variation, such as found by Rickard 
& Blitz (1985), could be attributed to variations in the ^"^CO abundance or excitation, 
rather than a spatial variation in X. Though separating these variations from other effects 
requires observations of other fines (e.g., ^^CO J = 3— >2) and isotopic variants (e.g., C^^O), 
a gradient in TZ would indicate an important variation in the physical or chemical properties 
of the molecular medium of a galaxy. 

In this paper, we present the results of a survey of the J = 1 — > fines of ^^CO and 
CO along the major axes of 17 galaxies (see Table 1). Unfike previous surveys, systematic 
uncertainties between intensities at different positions are reduced by using a receiver array. 
The array also permits efficient fully-sampled mapping of many positions in a galaxy. 



2. Observations and Data Reduction 

The data were obtained between 1998 December and 2000 May at the Five College Radio 
Astronomy Observatory (FCRAO) 14 m telescope. We observed the CO and ^^CO J=l— > 
lines (115.27 and 110.20 GHz) with the MMIC-based, 4 x 4-element array, SEQUOIA 
(Erickson et al. 1999). The beam FWHM at 115 GHz was 45", and the main beam 
efficiency was ~ 0.5. The pixel spacing on the sky was 1'476, or roughly two full beam 
widths. System temperatures were 130-400 K at 110 GHz, and 200-1200 K at 115 GHz. 
We used the facility broad-band filtcrbank, comprised of 64 channels of 5 MHz width, which 
yielded a velocity resolution of 13 kms^^ at 115 GHz. 

The observed galaxies are the brightest and most extended from the FCRAO Extra- 
galactic CO Survey of Young et al. (1995) that are primarily within the R.A. range of 
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9-16 hours (Table 1). The array was ahgned with the major axis of each galaxy, and moved 
from the edge of the CO disk across its full extent, with measurements taken at half-pixel 
(roughly half-beam) spacing. Thus each map position was observed with all four receivers of 
a row. Fully sampled maps were not made along the minor axes. Rather, off-axis data were 
obtained from the rows of the array parallel to the major axis, offset by ±1^476 and ±2.'952. 

Except for NGC 253, which has a low declination, we observed above an elevation of 
30° to avoid large gain variations. We also observed below an elevation of 75° to avoid 
any possible tracking errors. The pointing and focus were checked regularly on SiO maser 
sources and Jupiter, and calibration was done every several minutes using the chopper wheel 
method. A linear baseline was removed from all spectra. Those requiring higher order fits 
were discarded. Spectra with r.m.s. noise levels significantly different from those expected 
given the system temperature were also discarded. These data typically resulted from poor 
calibrations or non-linear baselines caused by quickly varying conditions. The filtered data 
were then coadded, weighted by the noise (1/a^). Therefore the largest calibration offsets, 
caused by low elevation and poor weather, contributed less to the overall results. Relative 
gain variations between pixels were well below 10%, and since each map position was observed 
by all four pixels in a row, this uncertainty was further reduced. 

To assess systematic uncertainties, the line integrated intensities were monitored for 
day-to-day variations (mostly at the central position). The weighted average CO and ^^CO 
intensities on any given day varied from the total season's weighted average by 5-25% and 15- 
50%, respectively. Any data with obvious ofi^scts in pointing (evident as markedly difi^crent 
spectral line shapes), or calibration (anomalous line and noise temperatures) were discarded. 
In sum, the relative spatial variations in the CO and ^^CO intensities, and their ratio, are 
very well determined for each galaxy. However, the systematic uncertainties in TZ, which are 
important for comparing the intensity ratios between galaxies, are roughly 15-60% (Table 2). 
Most of these systematic variations are likely due to slight pointing offsets along the minor 
axis, which are difficult to detect in low signal-to-noise spectra, and do not strongly affect the 
line shapes. Periodic observations of Jupiter indicate that the uncertainties due to variations 
in calibration contribute no more than 15-20% to the systematic uncertainty in line intensity. 



3. Results 

Emission from both CO and ^^CO is detected at nearly every mapped position along the 
major axes of these galaxies. Off-axis emission > lf5 away from the major axis, is detected 
only in the less inclined galaxies IC 342, M51 and NGC 6946. Due to the low signal-to-noise 
ratio of the ^^CO spectra off-axis, we do not include those data in the analysis. Figure 1 
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shows the spectra along the major axes of the survey galaxies. The similarity of the CO 
and ^^CO line shapes with position demonstrates the good relative pointing (better than 5") 
between the observations. The line integrated intensities (J T\dv) and CO/^^CO intensity 
ratios {TZ) as functions of position along the major axis are shown in Figure 2. 

The unweighted average of the CO/^^CO intensity ratio in the beam centered on the 
nucleus is 11.5 ± 0.3 ± 1.4, given the statistical and systematic uncertainties, respectively 
(Table 2). Statistical uncertainties are based on the r.m.s. noise, line width and channel 
width, and systematic uncertainties are estimated as described at the end of § 2. The 
standard deviation in TZ from one galaxy center to another is 5.4. The weighted average 
central ratios are 9.0 ±0.2 and 7.7±0.6 based on the statistical and systematic uncertainties, 
respectively. The central ratios for M51 and M82 deviate significantly from these values. The 
high ratio of 27 ± 8 for M82 is consistent with previous studies (e.g., Stark & Carlson 1982; 
Young & Scoville 1982; Loiseau et al. 1988; Kikumoto et al. 1998). However, the low ratio 
of 5.4±1.0 for M51 is nearly a factor of 2 below some previous observations (Young & Sanders 
1986; Garcia-Burillo et al. 1993). The individual CO and ^^CO data for M51 were taken over 
several widely separated days, and cover a large range in elevation, weather conditions, and 
times of day. Despite the varied conditions, the data sets agree well with one another (the 
estimated systematic uncertainties are 13 and 14% for CO and ^'^CO, respectively, yielding a 
19% uncertainty in TZ). The relative pointing was consistent, judging by the line shapes and 
radial intensity profiles. Therefore, we do not correct our values for M51 for any possible 
calibration offset, though we will draw no strong conclusions from the low ratio. We point 
out that our observing technique still results in a very low relative uncertainty in the spatial 
variation in TZ within the disk of M51. 

Figure 2 indicates that nearly all the galaxies show some significant variation in TZ with 
position. The most prevalent trend is for TZ to drop with galactocentric distance. The 
gradient in TZ is more apparent and significant after the line intensities are binned to make 
the uncertainty in TZ more uniform (Figure 3). Here non-detections have been excluded. 
The calculated uncertainties in these average ratios are multiplied by ^/2 to account for 
the interdependence of the Nyquist sampled data. This correction is conservative and may 
overestimate the error by about 20% for small bin sizes. Table 2 and Figure 4 compare the 
weighted average values of TZ inside and outside of various radii. Figure 4 shows that for 
many galaxies, TZ drops by as much as a factor of 2 between the central and outer disk. On 
average, TZ within a radius of 2 kpc is 40% higher than the disk value. No correlations are 
found between distance and TZ nor its gradient (Figure 4). 

Several galaxies show a significant change in the gradient in TZ with the size of the 
central bin (Figure 5). Note that typically the change in TZ does not depend strongly on the 
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central bin size, so the adopted distance to a galaxy is not very important. For those with 
significant changes, the bin size simply scales with distance, and intermediate values of TZ 
still lie on the trends shown. For NGC 253 and NGC 2146, the disk ratio decreases with 
the size of the central bin. The CO intensity distribution in NGC 253 is strongly centrally 
peaked inside 2 kpc. For NGC 2146, several outer points are not detected in ^^CO, so its 
distribution in the disk appears flatter than it may be (Figure 2), and the change in TZ is 
somewhat biased by this sensitivity hmit. For NGC 3556, NGC 4631 and NGC 6946, the 
central ratio decreases with the size of the central bin, and the gradient in TZ also decreases. 
NGC 3556 and NGC 4631 have flatter CO and ^^CO distributions with peaks offset from 
the center by 1.5-2 kpc that are more pronounced in ^^CO emission. The CO distribution 
for NGC 6946 peaks more sharply than the ^^CO at the nucleus. The central ratio increases 
with bin size for IC 342, NGC 3593 and NGC 3627. Their ^^CO distributions are more 
strongly centrally peaked within 1 kpc than the CO. M82 is unusual in that TZ in the disk 
increases as the central bin increases from 0.7 to 1.5 kpc radius, but then decreases for a 
central bin of 2 kpc radius. 



4. Discussion 

4.1. Possible Causes of Veiriations in the CO/^^CO Intensity Ratio 

Various physical mechanisms could cause a gradient in the observed CO/^^CO intensity 
ratio. The most obvious are changes in the relative CO and ^^CO abundances, and beam- 
averaged optical depths. A gradient in the X-factor, which may depend on the gas density 
and/or temperature, could also be responsible. For the following discussion, we assume, 
based on the centrally peaked proflles of integrated intensity (Figure 2), that the total beam- 
averaged H2 column density increases towards the nucleus of each survey galaxy. Also, we 
assume that TZ generally decreases with galactocentric distance. (It increases with radius 
only for NGC 3628, NGC 5055 and M51, and their individual ratios at large radii are still 
low.) 



4.1.1. Abundances 

There is a clear positive gradient in the ^^C/^'^C abundance ratio in the Milky Way (e.g., 
Langer & Penzias 1990, and references therein). It rises linearly from a value near 30 within 
4 kpc of the Galactic center, to ~ 70 at 10 kpc, and is most likely due to stellar processing 
of matter during the lifetime of the Galaxy; older populations have enhanced the ^^C in the 
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bulge. However, this gradient is in the opposite sense of the gradient in TZ seen here and in 
the Milky Way. For example, large-scale maps of the Milky Way suggest that TZ varies from 
10 ± 3 in the central 2 kpc to 6 ± 1 from 2 to 8.5 kpc (A. Luna, private communication), 
which is very similar to the results of this survey. Though TZ appears to be roughly constant 
at ~ 5 throughout most of the Galaxy (e.g., Solomon, Sanders & Scoville 1979; Polk et al. 
1988; Oka et al. 1998), large-scale observations may yield higher values of 7Z due to the 
inclusion of CO emission from very diffuse gas not detectable in ^^CO emission (Polk et al. 
1988; Lee, Snell & Dickman 1990). Also, the generally low CO/^^CO intensity ratios seen 
in galaxies indicate that at least the CO emission is optically thick, and therefore 7?. is a 
poor tracer of the ^^C/^^C abundance ratio. 

Another argument against stellar processing affecting TZ is that a correlation should 
exist between TZ and the C^^O/C^^O intensity ratio assuming that the stellar populations in 
the starburst nuclei of these galaxies contain relatively few low to intermediate mass stars 
(those responsible for much of the ^^CO generation) (Sage, Henkel & Mauersberger 1991). 
However, such a correlation is not observed. In particular, the C^^O/C^''0 intensity ratio 
for M82 of ~ 8 is similar to that of NGC 253 and IC 342, whereas the central CO/^^CO 
intensity ratio of M82 is twice those of the other galaxies. Casoli, Dupraz, & Combes (1992) 
also suggest that the CO abundance could increase by a factor of two due to an enhanced 
production of CO by massive stars. However, given the high optical depth of CO, this change 
will not alter TZ. Further, Taniguchi et al. (1999) attribute the high ratios seen in some 
galactic nuclei to a deficit of ^^CO, not an enhancement of CO. 

It has been suggested that the CO isotopic abundances may be altered by mechanical 
means as well. The high central CO/^^CO intensity ratios (> 20) seen in the nuclei of some 
mergers and ultraluminous starburst galaxies (Casoli et al. 1992a; Aalto et al. 1991; Henkel 
& Mauersberger 1993; Taniguchi et al. 1999), are thought to be due to either transport 
of less-processed (^^CO-poor) disk gas to the nucleus via an interaction, or destruction of 
dense (^^CO-rich) clouds by nuclear superwinds. If gas transport occurs faster than ^^CO 
production in the nucleus, any gradient should be erased. Models of galactic bars show that 
abundance gradients are indeed flattened after roughly 10^ yr (Friedli, Benz & Kennicutt 
1994), and that mass inflow time scales for slightly perturbed spirals are < 10^ yr, and shorter 
for stronger interactions (Jog & Das 1992, and references therein). That the gradient in TZ for 
the Milky Way opposes the ^^C/^^C abundance gradient also implies that altering abundances 
seems to have little effect on the CO and ^^CO emission. Further, NGC 3628 and M51 are 
two of the few galaxies in this survey with positive gradients in TZ, yet they are clearly 
in interactions. Therefore, unless a merger injects ^^CO-poor gas directly into a nucleus, 
transport due to bars or interactions seems unlikely to produce this gradient. The observed 
gradient in TZ instead favors the superwind argument since TZ and the star formation rate in 
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these galaxies are generally higher in the nucleus than in the disk. However, the volume and 
column densities of the gas toward these galactic nuclei are typically inferred to be large (e.g., 
Paglione, Jackson, & Ishizuki 1997). Therefore, if a superwind causes the higher CO/^^CO 
intensity ratio in the nucleus, it must force the clouds into many very small, dense clumps. 
This possibility is explored in the following section. 

Two chemical processes could be responsible for the gradient in TZ: isotope-selective 
photodissociation (ISPD) and chemical fractionation. ISPD can lower the ^'^CO abundance 
relative to CO in the nucleus since the rarer isotope is less shielded from destructive radiation. 
ISPD requires a strong, localized source of ultraviolet radiation, and dumpiness to maximize 
the surface area of the clouds while maintaining the high observed column densities. The 
nuclei of nearly all the survey galaxies exhibit high far-infrared (FIR) luminosities indicative 
of massive star formation, so sources of dissociating photons exist. ISPD would produce 
CO/C^^O intensity ratios even higher than TZ since the C^^O optical depth is lower than 
that of ^^CO. Nine galaxies have CO, ^^CO and C^^O emission observed at roughly the same 
resolution as this work (Casoli et al. 1992b; Aalto ct al. 1991; Sage et al. 1991), and there is 
a strong correlation between TZ and the CO/C^^O intensity ratio (Figure 6). However, there 
is no corresponding positive correlation between TZ and ^^CO/C^^O. In addition, ^^CO/C^^O 
stays the same or even decreases when viewed with higher resolution in NGC 253 and M82 
(Harrison et al. 1999; Wild et al. 1992) while it should presumably be highest near the 
central ionizing sources. In fact in M82 at 13" resolution, the highest ^^CO/C^^O ratio is 
found toward the weaker, eastern miUimeter continuum source (Wild et al. 1992; Carlstrom 
& Kronberg 1991). 

The C 11/ CO intensity ratio should be sensitive to the ultraviolet field strength and 
dumpiness of the gas (Stacey et al. 1991), and thus may also indicate the importance of 
ISPD. Figure 7 shows the central value of TZ for the survey galaxies plotted against this 
ratio. The linear correlation coefficient is less than 0.2, indicating a very low probability of 
correlation. 

Fractionation in the disk could lower the CO column density relative to that of ^^CO by 
factors of a few, but would have little effect on TZ because of the high optical depth of the CO 
line. Also, this mechanism requires lower temperatures (< 30 K), and CO optical depths 
< 10 (e.g., Keene et al. 1998). The CO line strengths indicate substantial disk column 
densities (> 10^° cm^^), but if the gas is sufficiently clumpy, fractionation on individual 
clump surfaces may be possible. Whether the cloud temperatures in the disk are low enough 
to support large-scale fractionation is more difficult to determine. The temperatures of the 
central regions of all the survey galaxies are inferred to be high from multi-line observations 
of CO (Wall et al. 1993; Wild et al. 1992), C 11 spectra (Stacey et al. 1991), and IRAS 
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colors (e.g., Young et al. 1989), but little data for galactic disks exist. Studies of NGC 6946 
and M51 indicate that higher gas temperatures in the disk arc limited to the spiral arms 
and various hot spots (Engargiola 1991; Fitt et al. 1992; Madden et al. 1993; Tuffs et al. 
1996). Little CO emission is seen outside of these warmer regions, where the gas might be 
cool enough to support fractionation. 

To summarize, we eliminate nearly all arguments in favor of abundance variations as 
the cause of a gradient in TZ. ISPD and CO enhancements due to gas transport or stellar 
nucleosynthesis are unlikely to be responsible for the observed gradients in TZ due to the high 
optical depth of CO, and the lack of expected gradients or correlations with other line ratios. 
Also the gradient in TZ in the Milky Way is in the opposite sense of the ^^C/^^C abundance 
gradient. It is difficult to test whether fractionation in the disk is important, though most 
of the CO emission in galactic disks is limited to relatively warm regions where it would not 
occur. Nuclear superwinds could effectively lower the beam-averaged ^^CO column density 
by creating very small clouds, which would result in higher values of TZ in starburst nuclei. 
This prediction is tested in the following section. 

4- 1-2. Cloud Properties 

The CO/^^CO intensity ratio should be a straightforward measure of the ^^CO optical 
depth. For example in LTE, given optically thick CO emission and optically thin ^^CO 
emission, and that their emission comes from the same volume of gas, 

7(CO)/7C=^CO) l/rC^CO). (1) 

In this case, the mean optical depth of the gas in these galaxies generally decreases toward the 
nucleus. In LTE, r oc {N/Av)/T^, for high T^. Therefore, to produce the observed gradient 
in TZ, the clouds in galactic nuclei arc probably warmer, have lower column densities, and/or 
higher internal velocity dispersions. 

To determine which cloud parameters dominate the variation in TZ, and to include vol- 
ume density as another important cloud parameter, we perform single-component model 
calculations of non-LTE CO excitation. We assume that the emission originates in unre- 
solved, homogeneous, spherical clouds, and include a photon escape probability function to 
account for the radiative excitation of optically thick hnes (Stutzki & Winnewisser 1985). 
The emission from the first 11 levels of CO is modeled. We use the coUision rates of Flower & 
Launay (1985). The excitation of ^^CO is assumed to be identical to that of CO (cf., Wilson, 
Langer, & Goldsmith 1981), and the ^^CO line temperature is estimated by decreasing the 
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column density of CO by the [COj/f^^CO] abundance ratio. Values of [CO] /[^^CO] =30-80 
are considered. 

Figure 8 shows the expected CO/^^CO intensity ratio from a cloud as a function of 
H2 volume density (nu^), CO column density per velocity interval (N/Av), and kinetic 
temperature (T)fc), given [CO]/[^^CO]=60. As expected, the ratio is a strong function of 
cloud column density and temperature, especially if the CO emission is thermalized. For 
example, a change in from 10 to 30 K results in roughly a factor 3 variation in TZ. For 
low densities, the optical depth of ^^CO is non-neglible, and the dependence on temperature 
diminishes (though a higher optical depth of the ^"''CO line does not necessarily imply 7^ ~ 1 
when LTE does not apply). Its dependence on density is weaker, especially for < 50 K, 
and reflects the difference between subthermally excited and thermalized (near LTE) CO 
emission. Below nH2~ 10^ cm~^, TZ may vary by factors of a few with density, depending on 
Tk and N/Av. 

To test how variations in N/Av may affect the gradient in TZ, we estimate the ratio 
of the column density within and outside of a 2 kpc radius in each galaxy. The column 
density per velocity interval for individual clouds, Nd Avc would have to decrease toward 
the nucleus to cause an increase in TZ. We also assume a proportionality between the CO 
and ^^CO intensities and the beam-averaged column density, N, such that 

(Arc/AVc)R>2kpc (^-^/ Av)R>2kpc ^ 0ii<2kpc ^ ' 

where (f) is the beam filling factor, and we assume constant CO and ^^CO abundances. The 
velocity width is estimated from the full line width. I / Av rises toward the nuclei of each 
survey galaxy by factors of 2-30 using CO, and 1.3-7 using ^^CO. For individual clouds, 
Nc/Avc could still decrease towards the nucleus if the beam filling factor rises faster than 
XI/ Av. In other words, if crowding in the nucleus is significant, then a higher value of TZ is 
still possible if Nc/Avc of individual clouds there is lower than in the disk. Such increases 
in (f) appear reasonable given the increase in line intensities and their likely dependence 
on varying excitation conditions from disk to nucleus. Therefore, we cannot eliminate the 
possibility of small nuclear clouds, such as would be produced by starburst winds, though 
the amount of variation in Nc/ Avc depends mostly on the internal turbulence of the clouds 
and the change in excitation conditions, that is, the variation in volume and column density 
and kinetic temperature. 

An increase in the gas volume density toward the nucleus can contribute to the gradient 
in TZ. For certain conditions {Tk > 10 K, nH2~ 10^-10^ cm~^), an increase in density of an 
order of magnitude can double TZ. The centers of galaxies are known to have enhanced gas 
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densities (Wall et al. 1993; Heifer & Blitz 1993; Paglione et al. 1997). Figure 7 shows 

the central value of TZ for the survey galaxies plotted against the (HCN J = 3^2)/ (HCN 
J = 1 — > 0) intensity ratio, which is sensitive to the gas density (Paglione et al. 1997). A 
positive trend is apparent, though there is no significant correlation (the linear correlation 
coefficient is 0.2). 

The correlation between TZ and IRAS color temperature (the ratio of global fluxes at 
60 and 100/im) is well established (Young & Sanders 1986; Aalto et al. 1991; Sage & Isbell 
1991; Taniguchi et al. 1999). There is a good correlation between TZ and color temperature 
for the galaxies presented here as well. According to the homogeneous cloud model, TZ 
varies with the gas kinetic temperature also, especially for > 10^ cm~^. Generally 7Z 
oc T^-^ for high densities, and TZ oc T^-^ for riHa < 10^ cm~^. Such a strong dependence 
on temperature, especially in combination with the expected increase in density towards 
the nucleus, suggests that most if not all of the variation seen in TZ is due to the high gas 
temperatures in these (mostly starburst) galactic nuclei. Figure 7 shows the central value 
of TZ for the survey galaxies plotted against the (CO J = 3 — > 2)/(C0 J = 1 — 0) intensity 
ratio, which is generally sensitive to T^. A linear correlation coefficient of 0.5 indicates a 
significant correlation between these ratios. 

A two-component model can also explain a rise in TZ at the nucleus without gradients 
in the bulk physical conditions of the clouds. We test a model where the cloud properties 
do not change, but warm and diffuse gas (nna = 100 cm~^, = 30 K) is added to the 
nucleus. Given the observed CO/^^CO ratios in galactic disks, a diffuse component with 
twice the beam filling factor of the dense component provides enhancements in TZ as high 
as 90% for very low temperatures (T^ < 20 K). Enhancements of 10-50%, for Tk = 20-75 
K and riu^ ~ 10^-10^cm~^, are more likely and still match the observed gradients in TZ. 
Decrements in TZ of up to 20-40% can also be produced by including the diffuse component, 
but require that the dense component be very dense and warm (riHa > 10^ cm~^, > 50 
K), which are unlikely conditions in galactic disks. 

To summarize, the observed gradient in TZ can be successfully explained by a gradient in 
the gas kinetic temperature. Increased volume densities in the nucleus may also contribute to 
raising TZ there. The column densities of individual clouds may decrease enough to cause the 
observed gradient in TZ, though this variation may be secondary to changes in temperature. 
This result is expected since in the LTE limit, the optical depth of ^"^CO, which is inversely 
proportional to TZ, depends more strongly on T^. Therefore a modest rise in towards the 
nucleus (e.g., from 10 K to 30 K, see Figure 8) can easily account for the gradient in TZ. 
Finally, a two-component model with warm and diffuse gas in the nucleus can reproduce the 
observed enhancements in TZ for a reasonable range of temperatures and densities. However, 



- 12 - 



the observed decrements in TZ are more difficult to explain with this model. 

4.1.3. X- Factor 

A positive gradient in the X-factor has been suggested for the Galaxy (Digel et al. 
1990; Sodroski et al. 1995; Dahmen et al. 1998). This result would follow naturally 
from the gradient in TZ found in this survey and the Milky Way. That is, the H2 mass 
estimated from a direct proportionality with CO intensity would overestimate the true mass 
in galactic centers, which is presumably traced less ambiguously by the more optically thin 
^^CO emission. Despite the probable Galactic gradient in the X-factor, the apparently strong 
dependence of TZ on cloud properties (and possibly chemical abundances) indicates that many 
things can contribute to varying the proportionality between Ico and H2 column density. 
Using the homogeneous cloud model from the previous section, as expected, TZ increases as 
X decreases given a constant CO abundance. Therefore, we suggest that a variable X-factor 
results from the same processes that affect TZ. Determining unambiguous variations in X 
requires a better understanding of the many physical processes that can affect TZ. 

We estimate the H2 column densities inside and outside a 2 kpc radius in the survey 
galaxies using various means: from the CO intensity and the standard, constant X-factor 
of 1.6x10^° cm~^ (K kms"^)"^; from the ^^CO intensity, assuming optically thin emission 
and LTE; and from the homogeneous cloud model, which includes non-LTE excitation and 
radiative trapping. 

The molecular masses inside a 2 kpc radius, derived from the standard X-factor, are 
hsted in Table 3. For most of the galaxies, the mass derived from CO is less than the 
dynamical mass inside 2 kpc. However for IC 342, M82 and NGC 6946, it is comparable to 
the dynamical mass. Unhke M82, IC 342 and NGC 6946 have typical CO/^^CO intensity 
ratios. 

The H2 column density can be estimated from the ^^CO intensity assuming optically 
thin emission and LTE using 

A^H2 =4.2xl0^9Tfee2-6^5/^^7(^=^CO), (3) 

where the ^'^CO abundance is 8xl0~^/60 (Frcrking, Langcr, & Wilson 1982). Masses derived 
from equation 3 are listed in Table 3 for = 10 K and 50 K. For most of the survey galaxies, 
these masses are similar to those derived from the standard X-factor and CO, which implies 
that the kinetic temperature is between 10 and 50 K. In fact, the gas kinetic temperature 
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can be estimated by setting the column densities derived from both CO and ^^CO. The 
temperatures inside 2 kpc (Table 3) are 20-60 K, and only somewhat lower in the disk. 
The higher temperatures are quite high compared to those derived from radiative transfer 
analyses (Wall et al. 1993), and seem unlikely over such large areas, especially in galactic 
disks. That this temperature seems high indicates that the masses derived from ^^CO are 
generally lower than those derived from CO. 

We estimate column densities with the homogeneous cloud model assuming = 300 
cm^'^ and = 10 K in the disk, and = 10'^ cm^'^ and = 30-50 K within a 2 kpc 
radius. Wc use the CO/^'^CO intensity ratios given in Table 2. Assuming LTE and optically 
thin ^^CO emission, the kinetic temperature is overestimated on average by a factor of 3-4 in 
the disk, and underestimated by 2.3-3 times in the nucleus. The average conversion factor in 
the disk is (l-3)xl0^°cm~^ (K kms~^)~^, which is consistent with the value in the disk of the 
Milky Way. For the nuclei we find X — (0.4-1.4) x 10^° cm~^ (K kms~^)"\ which is factors of 
2-5 times lower than in the disk, and 1.2-3.3 times lower than the standard Galactic value. 
Using the standard Galactic X-factor in galactic nuclei could therefore overestimate the true 
mass by factors of a few. Note that the model indicates that the X-factor is always smaller 
in galactic nuclei, regardless of the gradient in TZ. 

5. Conclusions 

We observed the CO and ^^CO J = 1 ^0 emission along the major axes of 17 nearby 
galaxies. On average, the CO/^^CO intensity ratio, TZ, inside a 2 kpc radius is roughly 
30% higher than in the disk. This ratio is sensitive to variations in temperature and column 
density, as well as fractionation and isotope-selective photodissociation. Winds and gas inflow 
caused by mergers have also been suggested to vary TZ, along with abundance variations due 
to stellar processing. We eliminate most mechanisms except for fractionation and variations 
in kinetic temperature, cloud column density and H2 volume density. The gradient is most 
likely caused by the higher temperatures (and perhaps densities) typical of the central regions 
of starburst galaxies. Small nuclear clouds, perhaps caused by starburst superwinds, may 
also contribute to the variation in TZ. We estimate that the X-factor (X — N-a^/Ico) 
decreases toward galactic nuclei, as seen in the Milky Way, and its variation results from the 
same physical processes that affect TZ. A modest increase in global gas temperatures easily 
accounts for the observed variation in TZ, as well as a decrease in X of factors of 2-5, from 
disk to nucleus. 
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Table 1. Properties of Survey Galaxies and Maps. 



Galaxy 


R.A. (2000) 


Decl. (2000) 




D 


Map P. A. 




Map Size'^ 




h m s 


o / // 


(kms~^) 


(Mpc) 


(degrees) 


(kpc) 


(kpc) 


NGC 253 


00 47 35.1 


-25 17 20 


230 


3 


54 


0.7 


9.3 


NGC 1068 


02 42 40.7 


-00 00 48 


1150 


UA 


90 


3.2 


20.1 


IC 342 


03 46 49.7 


68 05 45 


30 


4 





0.9 


15.9 


UGC 2855 


03 48 22.6 


70 07 57 


1200 


20 


92 


4.5 


30.1 


NGC 2146 


06 18 37.6 


78 21 19 


920 


14 


128 


3.1 


13.5 


M82 


09 55 54.0 


69 40 57 


250 


3.25 


65 


0.72 


8.4 


NGC 3079 


10 01 58.2 


55 40 42 


1150 


20 


164 


4.5 


27.9 


NGC 3184 


10 18 17.2 


41 25 26 


590 


13 


135 


2.9 


23.7 


NGC 3556 


11 11 31.8 


55 40 14 


700 


10 


80 


2.2 


15.0 


NGC 3593 


11 14 36.0 


12 49 06 


600 


9 


90 


2.0 


12.6 


NGC 3627 


11 20 14.4 


12 59 42 


740 


6.7 





1.5 


12.2 


NGC 3628 


11 20 16.2 


13 35 22 


850 


6.7 


103 


1.5 


15.1 


NGC 4527 


12 34 08.8 


02 39 13 


1770 


15 


67 


3.3 


20.9 


NGC 4631 


12 42 07.6 


32 32 28 


630 


9 


88 


2.0 


13.5 


NGC 5055 


13 15 49.2 


42 02 06 


500 


7 


105 


1.6 


15.8 


M51 


13 29 53.2 


47 11 48 


450 


9.6 





2.1 


21.6 


NGC 6946 


20 34 51.8 


60 09 15 


50 


10 


45 


2.2 


22.5 



^Central position 

''Linear extent of 46" at distance D 

■^FuU extent of map along major axis, sampled every 22'.'14 with a 46" beam (FWHM). 
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Table 2. Variation of CO/^^cQ Intensity Ratio. 



Galaxy 




I{CO)/lC^COy Ratio'' 




R < 23" 


R > 23" i? < 2 kpc i? > 2 kpc 



NGC 253 


11.9 


± 


0.5 


± 


3.4 


8.8 


± 


0.5 


12.4 


± 


0.5 


6.9 


± 


0.8 


1.8 


± 


0.2 


NGC 1068 


10.7 


± 


0.6 


± 


1.9 


11.7 


± 


1.2 


10.4 


± 


0.4 


11.6 


± 


1.1 


0.9 


± 


0.1 


IC 342 


8.4 


± 


0.5 


± 


1.6 


7.1 


± 


0.6 


10.8 


± 


0.8 


6.4 


± 


0.6 


1.7 


± 


0.2 


UGC 2855 


11.7 


± 


1.7 


± 


2.3 


9.4 


± 


1.7 


11.7 


± 


1.7 


9.4 


± 


1.7 


1.2 


± 


0.3 


NGC 2146 


14.9 


± 


1.3 


± 


8.6 


28.1 


± 


6.9 


15.5 


± 


1.1 


22.8 


± 


4.5 


0.7 


± 


0.1 


M82 


27.3 


± 


1.0 


± 


7.8 


12.9 


± 


0.5 


20.9 


± 


0.7 


12.1 


± 


6.9 


1.7 


± 


1.0 


NGC 3079 


15.3 


± 


1.8 


± 


5.1 


10.1 


± 


2.1 


15.3 


± 


1.8 


10.1 


± 


2.1 


1.5 


± 


0.4 


NGC 3184 


4.9 


± 


1.2 


± 


2.3 


4.4 


± 


0.8 


5.6 


± 


1.1 


4.0 


± 


0.6 


1.4 


± 


0.3 


NGC 3556 


8.7 


± 


1.4 


± 


5.6 


6.8 


± 


0.6 


7.0 


± 


0.8 


6.9 


± 


0.7 


1.0 


± 


0.2 


NGC 3593 


10.6 


± 


2.2 


± 


10 


11.4 


± 


3.2 


16.8 


± 


4.6 


5.7 


± 


2.3 


2.9 


± 


1.4 


NGC 3627 


11.2 


± 


1.2 


± 


5.7 


17.4 


± 


2.1 


15.5 


± 


1.6 


17.8 


± 


4.1 


0.9 


± 


0.2 


NGC 3628 


9.2 


± 


0.4 


± 


2.3 


10.3 


± 


0.6 


9.7 


± 


0.4 


11.2 


± 


1.1 


0.9 


± 


0.1 


NGC 4527 


6.1 


± 


0.5 


d= 


1.4 


6.9 


± 


0.8 


6.7 


± 


0.9 


6.7 


± 


0.6 


1.0 


± 


0.2 


NGC 4631 


16.3 


± 


2.4 


± 


9.9 


9.1 


± 


0.7 


13.1 


± 


1.3 


8.4 


± 


0.9 


1.6 


± 


0.2 


NGC 5055 


5.7 


± 


0.6 


± 


3.1 


8.2 


± 


0.7 


7.1 


± 


0.6 


9.2 


± 


1.1 


0.8 


± 


0.1 


M51 


5.4 


± 


0.4 


± 


1.0 


7.5 


± 


0.6 


6.4 


± 


0.4 


8.2 


± 


0.9 


0.8 


± 


0.1 


NGC 6946 


17.0 


± 


1.4 


± 


6.6 


10.4 


± 


1.1 


12.7 


± 


1.0 


9.8 


± 


1.4 


1.3 




0.2 


Average 


11.5 


± 


0.3 


± 


1.4 


10.6 


± 


0.5 


11.6 


± 


0.4 


9.8 


± 


0.6 


1.3 


± 


0.1 



^Prom weighted averages of data within or beyond noted radii. For the cen- 
tral position, the estimated systematic uncertainty is listed after the statistical 
uncertainty 

^Intensity ratio inside 2 kpc divided by that outside of 2 kpc. Uncertainty is 
based on the statistical uncertainty. 



Table 3. Modeling Results". 



Galaxy 


Mh,(co) Mh,(13co) n 


(K) 


X 




Tfc = 10 K Tfc = 50 K i? < 2 kpc 


i? > 2 kpc 


i? < 2 kpc i? > 2 kpc 



NGC 253 


2.40 


(0.04) 


0.66 


(0 


03) 


2.7 


(0 


1) 


45 


(2) 


23 


(3) 


n r+0 4 


Z.O_Q 7 


NGC 1068 


2.16 


(0.32) 


0.75 


(0 


-1 r\\ 

10) 


3.0 


(0 


4) 


37 


(2) 


42 


(4) 


o-7io:2 


1 1 +0 4 


IC 342 


0.99 


(0.16) 


0.32 


(0 


03) 


1.3 


(0 


1) 


39 


(3) 


22 


(2) 


^ „+0 3 

o-6io;2 


„ o+O 9 


UGC 2855 


0.87 


(0.03) 


0.25 


(0 


04) 


1.0 


(0 


1) 


42 


(7) 


33 


(6) 


6+°-3 


1 6+°-5 

-'-•0_o.4 


NGC 2146 


1.47 


(0.29) 


0.34 


(0 


09) 


1.4 


(0 


4) 


56 


(4) 


84 


(17) 


4+^-° 




M82 


2.83 


(0.01) 


0.46 


(0 


01) 


1.9 


(0 


1) 


77 


(3) 


44 


(26) 


0.3+°;? 


-'-•-'--0.4 


NGC 3079 


1.87 


(0.02) 


0.42 


(0 


05) 


1.7 


(0 


2) 


56 


(7) 


36 


(8) 


o.4t°:^ 


i-4i^i 


NGC 3184 


0.16 


(0.01) 


0.09 


(0 


02) 


0.4 


(0 


1) 


19 


(4) 


13 


(2) 


1 5+2'^ 
-'-■'J-0.7 


4.9l?i 


NGC 3556 


0.34 


(0.01) 


0.17 


(0 


02) 


0.7 


(0 


1) 


24 


(3) 


23 


(3) 


I 1+5.0 


9 r+8.7 
^■^-1.3 


NGC 3593 


0.50 


(0.04) 


0.10 


(0 


03) 


0.4 


(0 


1) 


61 (17) 


19 


(9) 


> 0.1 


> 1.2 


NGC 3627 


0.81 


(0.01) 


0.18 


(0 


02) 


0.7 


(0 


1) 


57 


(6) 


65 


(15) 




U.O_Q 2 


NGC 3628 


1.47 


(0.02) 


0.52 


(0 


02) 


2.1 


(0 


1) 


34 


(2) 


40 


(4) 


o.7;gi 


-1 9+0.6 
-'-•^-0.4 


NGC 4527 


1.03 


(0.10) 


0.49 


(0 


12) 


2.0 


(0 


5) 


23 


(3) 


23 


(2) 


1 9+0.8 
-"^•^-0.4 


2 6+^-2 

Z.O_Q 7 


NGC 4631 


0.83 


(0.02) 


0.22 


(0 


02) 


0.9 


(0 


1) 


47 


(5) 


29 


(3) 




-'-•8_i.o 


NGC 5055 


0.57 


(0.02) 


0.27 


(0 


02) 


1.1 


(0 


1) 


25 


(2) 


32 


(4) 


1 1+2-^ 


-'-•"-0.8 


M51 


1.16 


(0.02) 


0.62 


(0 


04) 


2.5 


(0 


2) 


22 


(2) 


29 


(4) 


i.3;S;I 


1 q+0-7 

-■-•^-0.4 


NGC 6946 


1.79 


(0.03) 


0.48 


(0 


04) 


2.0 


(0 


1) 


46 


(4) 


35 


(5) 




1 5+^-5 



'^Uncertainty given in parentheses. 
Column 1: Galaxy 

Column 2: Mass inside 2 kpc from standard conversion factor and /(CO) in 10^ Mq 

Columns 3 and 4: Mass inside 2 kpc from ^^CO (eq. 3) for given in 10^ Mq 

Columns 5 and 6: from standard X-factor and column density estimated from ^^CO (eq. 3) 

Columns 7 and 8: X-factor from non- LTE calculations and 7^ (Table 2) in 10^° cm-2(K 
kms-^)-i 
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Fig. 1. — Spectra of CO (heavy lines) and ^^CO (light hnes) emission for the survey galaxies. 
Except for IC 342, all ^^CO spectra are smoothed to 27.2 kms~^ resolution for display. The 
temperature scale is for the ^^CO line. The CO line is divided by 10 for display. 

Fig. 2. — The CO and ^^CO integrated intensities (jT^dv), and their ratios, versus position 
along the major axis. Upper and lower limits (Sa) are denoted with downward- and upward- 
pointing triangles, respectively. Error bars are la statistical estimates based on the r.m.s. 
noise, line width and channel width. 

Fig. 3. — The CO/^^CO intensity ratio as a function of position. The data are folded about 
the nucleus and binned to make the uncertainty in TZ uniform along the major axis. Vertical 
error bars indicate the statistical uncertainties of the data (see text). Horizontal bars indicate 
the bin sizes. 

Fig. 4. — Comparison of the ratio of CO and ^^CO integrated intensities inside and outside 
of selected radii: 23" (beam HWHM), 1, 1.5 and 2 kpc. The hnes indicate slopes of 2/3, 1, 
3/2, 2 and 4 (dash-triple-dot, sohd, dashed, dash-dot and dotted, respectively). 

Fig. 5. — Same as Figure 4 indicating the variation in the gradient of TZ with the size of the 
central region. The ratios for M82 are divided by 10 for clarity. The filled symbol denotes 
the ratio for the central point without binning. The other connected points are for central 
regions of 1, 1.5 and 2 kpc radii, in sequence. 

Fig. 6. — Central CO/^'^CO intensity ratio plotted against the central CO/C^^O and 
"'^^CO/C-'^^O intensity ratios. All data arc at similar resolution, ~ 45". References - NGC 253, 
IC 342, M82: this work, Sage ct al. (1991); NGC 1808, NGC 2146, NGC 4826, Circinus, 
NGC 7552: this work, Aalto et al. (1991); NGC 3256: Casoh et al. (1992b). 

Fig. 7.— Central CO/^^CO intensity ratio plotted against C II/CO (Staccy et al. 1991; 
Carralet al. 1994), HCN J = 3^2/J=1^0 (Paglione et al. 1997), and CO J = 3^2/J = 
1 — s> (Mauersberger et al. 1999; Papadopoulos & Seaquist 1999; Wielebinski, Dumke & 
Nieten 1999). Error bars denote systematic uncertainties. 

Fig. 8. — Expected CO/^^CO intensity ratios as functions of kinetic temperature, CO column 
density per velocity interval, and H2 volume density. 
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